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A rocket pump inducer is usually an axial-flow pump runner hav- 
ing very high solidity and very low aspect ratio blades operating at 
very low coefficients. While this characteristic form is dictated by 
cavitation requirements, the flow is subjected to major effects of 
viscosity and turbulence in the flow field in the long and narrow pass- 
ages between the vanes. The investigations reported in this paper 
are concerned with the effects of viscosity and not with the effects of 
cavitation. 

The experimental investigations are carried out on a 3-foot-diam- 
eter model of a four-bladed inducer which is operated in air at a flow 
coefficient 0.065. The fluid properties are measured at the exit of the 
inducer using conventional and hot-wire probes. Important experi- 
mental results and the method of predicting the outlet tangential ve- 
locity and head rise are discussed in this paper. 


One of the major requirements of a turbopump used in rocket applica- 
tion is that it should be capable of running at very high speeds, so as to 
minimize the size and weight of the unit and to facilitate matching with a 
drive turbine. The conventional pumps cavitate at suction specific speeds 
(SS) in excess of 8000, thus limiting the highest possible speed. The need 
to increase the speed led to the development of a cavitation resistance 
inducer which is essentially an axial-flow pump with high solidity blades 
used in front of the main pump. Long and narrow helical passages provide 
a space and time for the collapse of cavitation bubbles and for the gradual 
addition of energy to the fluid. Inducers are also used in water jet pumps 
of marine application (ref. 2). These are light and feature cavitation 
resistance in excess of commercial applications. 


1 The work reported here was done under contracts NSG537 and NGL 39-009-007 
of the National Aeronautics and Space Administration, with technical management 
by W. R. Britsch of NASA Lewis Research Center. 
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The purpose of the inducer is to pressurize the flow sufficiently to enable 
the main pump to operate satisfactorily. The physical reasoning for the 
selection of such unconventional blade passages is explained by Acosta 
(ref. 1). These inducers have operated successfully at suction specific 
speeds in excess of 30 000; the test inducer described later is designed for 
SS = 50 000. The characteristic features of the inducers are: 

(1) High solidity blading (few long-chord blades instead of many 
short-chord blades as in axial-flow pumps and compressors, exposing a 
much larger area to frictional effects) 

(2) Very low flow coefficient and large stagger angles. 

While the characteristic form of the inducer is dictated by cavitation 
requirements, the flow is subjected to major effects of viscosity and 
turbulence of the fluid in the long and narrow passages between the 
inducer vanes. The investigations reported here are concerned with the 
effects of viscosity and not with the effects of cavitation. 

In axial-flow compressors or pumps, the secondary motions or the 
departure of the flow from design values are confined to fairly thin regions 
around the blades and near the end walls. This is not true in the case of 
inducers used in liquid rocket pumps. The blades of such inducers extend 
circumferentially over a major portion of a complete circle (fig. 1). Under 


Figure 1 . — Photograph of the inducer ( without shroud); 36-inch outside diameter, four 

blades. 
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these circumstances, it can bo estimated that the viscous effects and the 
consequent secondary motions are many times stronger than those en- 
countered in other types of axial-flow machinery, and may indeed domi- 
nate the entire flow field. An understanding of such secondary motions is 
essential for the prediction of the efficiency and general performance of the 
inducer. The radial distribution of flow properties and the appropriate 
method of predicting them provide a further step toward the under- 
standing of the flow through an inducer. This type of information is 
necessary for the analysis and design of the main pump. This is the 
objective of this paper. 

The Department of Aerospace Engineering at The Pennsylvania State 
University has undertaken a systematic investigation with a view to 
understanding the general flow behavior and to developing flow pre- 
diction methods for inducers. A 3-foot-diameter inducer operated in air at 
a flow coefficient of 0.065 was built for this purpose. The following in- 
vestigations, using the test rig shown in figure 2, have been completed 
so far: 

(1) Visualization and measurement of the flow inside the passages* 
and at the exit of a four-bladed inducer. The latter measurements and 
an approximate analysis arc reported in this paper. 



NASA In * 


Figure 2. — Test assembly. 


2 These are carried out using rotating probes and pressure transfer device described 
in reference 20. The measurements are reported in reference 21. 
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(2) Three-dimensional velocity and turbulence measurements at the 
exit of a three-bladed inducer using hot-wire probes. Cooper and Bosch’s 
computer program (ref. 9) is used to predict the three-dimensional inviscid 
effects (refs. 3 and 4). 

(3) Three-dimensional turbulent boundary investigation on a single 
helical blade of the type used in rocket pump inducers (ref. 5). 

These and future investigations to be undertaken will lead to the 
establishment of a theoretical model for the eventual analysis and design 
of turbomachinery flow dominated by secondary motions. 

It is the intention of this paper to discuss the exit flow characteristics, 
a correlation for the frictional losses for the inducers, and an approximate 
analysis for the predictions of the velocity and energy distribution at the 
exit of the inducer. It will be shown that, while adequate information is 
not available to develop an exact analysis, the predictions based on loss 
correlations provide a method of predicting the gross behavior of the 
inducer. 


DESIGN OF TEST INDUCER 

The design of the test inducer is based on the “mean streamline” 
method developed by Wislicenus (ref. 6). The detailed design is described 
in reference 7. The overall characteristics arc 1 listed in table I. 

After selecting a suction specific speed (SS) of 50 000 and head coeffi- 
cient \p T = 0.2, the following parameters were derived from the design 
chart relating the cavitation parameters of turbomachines (fig. 259 in 
reference 6) . 

Flow coefficient, <£ = 0.065 

Blade minimum pressure coefficient, Cp min = 0.0 1 

2 (jH sv /V z 2 

The runner is of free vortex design, producing constant head from root 
to tip. Using the relationship relating the stagnation head coefficient to 
change in tangential velocities, the inlet and outlet velocity diagrams are 
derived (fig. 3). The resulting runner profile is shown in figure 4. The 
design distribution of tangential and axial velocity components and of the 
stagnation and static head coefficients are shown in figures 11, 16, 18 and 
20, respectively. 

The design of blade sections is carried out in three steps: 

(1) The mean streamline is derived from the assumed blade pressure 
diagram and blockage distribution as shown in figure 3. The blade pres- 
sure distribution chosen is typical of the trailing edge loaded profile which 
is required for this type of turbomachine and has Cp min = 0.01. In order to 
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account for the thickness or blockage effect of the blades, the relative 
velocities are varied from Wi to W 2 along an arch-shaped curve above the 
nearly straight line variation that would apply to zero blade thickness. 
The arch-shaped curve represents an estimate of the blockage due to 
blade and boundary layer displacement thickness. The assumed blockage 
due to both effects is clearly marked in figure 3. The mean streamline for 
the hub and tip sections derived from the assumed pressure diagram and 
blockage is shown in figure 3. 

(2) The departure of the camber line from the mean streamline is 
derived from the empirical deviation values of A rj/L (nondimensionalized 
distance between the mean streamline and camber line measured normal 
to blade chord) derived by Wislicenus (ref. 6). It should be remarked here 
that Wislicenus’ correlation for A i\/L is the empirical counterpart of 
Ackeret’s theoretical method (ref. 8). The 1 validity of these empirical 
correlations for standard profiles has been established by Jakubowski 
(ref. 6) but its generalization to passages with nonstandard blade profiles 


Table I . — Inducer Characteristics 


Tip diameter. __ .. 


36.5 inches 

Hub /tip ratio at outlet 



0.5 

Hub/tip ratio at inlet . .. 


0.25 

Radial clearance 



0.0625 inches 

Suction specific speed SS (design) 

'. _ 

. 50 000 

Flow coefficient (4> = V,/Ui). . . 



0.065 

Blade chord. . . 

. Tip (ff = 1.0) 

_ 82.06 inches 


Mid span (ft = 0.75).. 

63.18 inches 


Hub (ft = 0.5) 

49.04 inches 

Solidity. . . ... 

. Tip (ft = 1.0) 

. 2.86 


Mid span (ft =0,75). . . 

2.01 


Hub (ft =0.5) 

. . 3.50 

Number of vanes - 


4 

Angular wrap (average). 

Lift coefficient of the blade based 


_ 290 c 

on mean velocity . ... 

. Tip 

0.0966 


Mid span .... 

_ 0.163 


Hub 

0.307 

Reynolds number based on tip radius.. 
Reynolds number based on relative 


6.60 X10 5 

velocity and chord at mid radius. 
Maximum deviation of camber line 


1.75X10 6 

from mean streamline) Arj/L) ma:t 

(ft = 0.5) 

_ 0.02 


(ft =0.75) 

0.01075 


(ft = 1.0) 

_ 0.00637 

Blade angles at inlet 

.. (ft =0.25) 

. _ 75°30' 


(ft =0.625) 

. _ 83°30' 


(ft = 1.0) 

86 T5' 


All other blade angles are given in figure 0. 
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Figure 3. — Design of blade profiles. 


SECTION 




Figure 4. — Axial view and circular projection of the inducer blade. 

and large stagger angle and solidities is yet to be established. The chord- 
wise variation of Ai j/L assumed for this design departs from the XACA 
65-series cascade data as required to obtain trailing edge loading (fig. 5). 
In figure 5, the notation (Ay/L) m w, 65 refers to the maximum deviation for 
a NACA 65-series cascade. Its value can be evaluated from the equation 



max, 65 


max 
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Figure 5. — Assumed deviation of the camber line from the mean streamline. 


where 


C L lift coefficient 


maximum deviation for unit lift coefficient 

Their values for various vane chord angles are plotted by Wislicenus 
(fig. 298, ref. 6). The value of (Anji/ L) max used for the inducer is 0.05 and 
the lift coefficients at various radial sections are given in table I. The 
camber lines so derived at hub and tip sections are shown in figure 3 and 
the design blade angles derived from this method are shown plotted in 
figure 6. 

(3) The blade thickness can be derived from the assumed blockage 
curve (fig. 3) and the equation 

A Vz_t 

~V7~s 



where 

AFz perturbation in local axial velocity V z due to blade blockage 
t blade thickness in tangential direction 
s blade spacing 

The blade profiles so derived at hub and tip sections are shown plotted 
in figure 3. 
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Figure 6. — Design blade and air angles. 



NON-DIMENSIONALIZED RADIUS, R 


The blade section at hub (72 = 0.5) , midradius {It = 0.75) and tip 
sections (/£=1.0) are designed in the manner explained above. The 
stacking of blade sections is carried out in much the same way as described 
by Wislicenus (ref. 6) and the final blade system so derived is shown in 
figure 4. In order to achieve reasonably constant variation in blade thick- 
ness from root to tip and smooth curvature in both chordwise and radial 
directions, it is sometimes necessary to change blade sections. This re- 
sulted in a slight departure from the free vortex design, especially at mid- 
radius, as shown in figures 20, 18, 16 and 11. 

It should be remarked here that the Penn State inducer differs appre- 
ciably from th(> conventional radial blades of constant thickness (usually 
called flat plate inducer) used in practice. From the point of view of 
manufacture the radial blades are superior, but the design based on the 
present method should have hydrodynamically superior performance. 

The blades for the inducer are of fiberglass construction, cast from 
suitable mold and templates. The final assembly of the test inducer, which 
is driven by a 5-lip variable-speed motor, is shown in figures 1 and 2. 
The research rotor is fitted with a pressure transfer unit which transmits 
pressures from a rotating blade to a stationary manometer. 

APPROXIMATE THEORETICAL ANALYSIS 

An attempt has been made to predict the flow properties at the exit of 
the inducer using a circumferentially averaged radial equilibrium equa- 
tion. The three-dimensional nature of the flow and energy losses due to 
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friction are fully recognized in deriving the circumferentially averaged 
radial equilibrium equation for the inducer. The inducer passages are 
assumed to have fully developed turbulent flow. Based on the information 
available in the current literature, assumptions are made for the circum- 
ferential variation of the relative tangential and radial velocities. The 
frictional losses accounted for in the theory are based on a friction loss 
coefficient (analogous to the universal friction loss coefficient derived by 
Blasius for pipes) derived from various inducer tests carried out at 
NASA Lewis Research Center, TRW Cleveland, and M.I.T. Gas Turbine 
Laboratory. This friction loss coefficient, which takes into account 
rotational effects, is valid for low flow coefficient (upstream absolute 
velocity/blade tip speed) or high rotation parameter (blade tip speed/ 
upstream absolute flow). The energy losses calculated from this newly 
derived friction loss coefficient are in close agreement with the energy loss 
measured in the Penn State inducer. The radial equilibrium equation 
derived in this section is valid for a noncavitating inducer. 

The various assumptions made in deriving the radial equilibrium 
equation for the inducer are given in the text in proper context. The 
absolute and relative tangential velocities predicted from this equation 
are in good agreement with the measured values (see “Experimental 
Results and Discussion”) . The velocity profile models used in this analysis 
are based on the three-dimensional boundary layer investigation carried 
out by Lakshminarayana et al. (ref. 5) for a single blade. The viscid effects 
art 1 taken into account in an approximate manner by using an empirically 
derived loss coefficient in the final radial equilibrium equation. The analy- 
sis considers only the perturbations caused by viscous effects. The inviscid 
turning effects are allowed for by taking design or the potential flow solu- 
tion as the boundary conditions for the solution of the radial equilibrium 
equation presented in this section. 

No attempt has been made to predict the axial velocity components, 
which are very much smaller than the tangential velocity components. 
This prediction should await more accurate information on the three- 
dimensional nature of the flow in the rotating passages. 


Equations of Motion 

The equations of motion governing the incompressible and steady flow 
inside the passages in a coordinate system rotating with an angular 
velocity (fig. 7) are given by reference 19, 

V/ = 2WX£2+WX(VXW)+F (1) 


V-W = 0 


( 2 ) 
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where 


7 = 


P 

-+ 

p 


W 2 

2 


Q 2 r 2 

IT 


(3) 


and F is the frictional force per unit mass. 

Since large velocity and pressure gradients exist in the case of an in- 
ducer, it is a formidable or impossible task to obtain an exact solution of 
equations (1) and (2) . Hence, a more practical approach to satisfying the 
radial equilibrium equation on a circumferential average basis is employed 
in this paper. This technique was first employed by Ruden (ref. 10) for a 
compressor and later used by Smith (ref. 11) to derive a radial equilibrium 
equation amenable to solution by what is known as the “stream line 
curvature approach.” 

The scalar component of equation (1) in the radial direction is 


— ^-We^-W 

dr rdd 


dW r , „ r dWe dWz 
— — {-We - , 

dz dr dr 


We 


- 2 OIT 9 + 7 1 T (4) 


The following assumptions are made in deriving a radial equilibrium 
equation applicable to inducers operating at low flow coefficients. 

(1) Frictional effects: The fluid friction affects the radial equilibrium 
equation through the terms F r and 7 in equation (4). The radial com- 
ponent of shear stress at the blade surface is small, as indicated by 
Lakshminarayana’s investigation (ref. 5) on a helical blade. It is recom- 



Figure 7. — Notations used for inducer flow analysis. 
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mended that the shear stress effect F r be neglected in the preliminary- 
analysis. In an inviscid flow, I is constant along a streamline. In viscid 
flow it can be approximated by 


I 


p ip nv aPq 

p ' 2 2 ’ p 


where A P 0 /p represents the stagnation pressure losses along a streamline. 
This term includes all losses associated with friction and turbulent mixing. 
In long, narrow inducer passages, this term is likely to be large. Thus, 


dj__ _d^ /APiA 
dr dr\p) 


( 5 ) 


(2) For inducers operating at low flow coefficients, 

W z «W e 

OWz dll e 
dr dr 


(6) 


(3) The flow passages arc assumed to have fully developed flow and 
the flow properties are symmetrical about the midpassage (fig. 7). The 
assumptions for boundary velocity profiles are based on the investigations 
of reference 5 and are given by equations (12) and (13). 

(4) The fluid is assumed to be guided smoothly through the channel 
(i.c., no flow deviation) . 

(5) The fluid is assumed to be single' phase, incompressible and steady, 

Assumptions (1) and (2) above would reduce equation (4) to 


w dlFr dWr 

rdd dz 


TI7 <9 IF, IIY 
-IF« - — I 

dr r 


d /APo 

-2QIF 9 = 

dr\ p 


( 7 ) 


1 2 3 4 5 

Prior to deriving a circumferential average of the above equation, an 
attempt is made to present a method of evaluating the stagnation pressure 
loss (in relative flow) from inducer loss correlations. 


Friction Loss Coefficient for a Rotating Channel with Large 
Rotation Parameter or Small Flow Coefficient 

Since relative velocities are zero at the solid boundaries of a rotating 
channel, one is tempted to contemplate the behavior of the boundary 
layer in a frame of reference fixed to the rotor. 

The fluid near the solid blade surfaces is not subjected to the same 
Coriolis or centrifugal forces as the main flow, and the result is that 
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externally impressed gradients produce secondary motions inside the 
blade boundary layer. The additional losses associated with these second- 
ary motions increase the frictional losses. Investigations by Spannhake 
(ref. 13) and Ludweig (ref. 14) indicate that the loss coefficient in such 
rotating channels may be several times that of an equivalent 3 stationary 
channel. Very little is known about the increase in friction losses due to 
rotation, especially the additional frictional stresses induced by Coriolis 
force. Systematic experimental investigations by Spannhake (ref. 13) and 
Ludweig (ref. 14) indicate that these frictional losses are strongly de- 
pendent on the rotation parameter. This parameter is analogous to the 
Rossby number used in connection with studying the effect of earth 
rotation on wind profile and is defined as the ratio of Coriolis to inertia 
force. 

Since very little analytical information is available for the type of 
rotating channel used in the inducer, a systematic attempt is made to 
correlate the measurements of various investigators to derive a friction 
factor for inducer channels for the range of rotation parameters (inverse 
of flow coefficient) used in practice. The inducer data collected are from 
various laboratories under different flow conditions and blade configura- 
tions. The frictional losses depend on Reynolds number, velocity, dimen- 
sions of the channel, rotation parameter, and aspect ratio of the blades. 
Since the available experimental data on inducers are meager and the 
tests involve a small range of Reynolds number, no systematic approach 
is used to find the Reynolds number dependency on losses. First, the 
Blasius friction factor (X) is determined from the equation 


APo_ X L _ 2 
p 2R N l/i d h 


( 8 ) 


where 

AZ-’o/p is the loss in stagnation pressure of the relative flow 
X is the Blasius friction factor (0.310 for a stationary channel) 

R N is the Reynolds number, Wdjv (fig. 7) 

dh is the hydraulic mean diameter (4 X area/wetted perimeter) 

L is the length of the channel 

For an elemental blade height dr of the inducer passage at any radius 
r (fig. 7) 


dh = 


4(2tt r/n) cos (3 dr 
2 dr 



COS d 


( 9 ) 


where (3=f{r) and n is the number of passages. 


3 Equivalent here refers to a fully developed stationary channel having a mean 
velocity W. 
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It is easy to prove that the nondimensionalized stagnation head loss 
coefficient in relative flow is equal to the difference between the Euler and 
measured head coefficient in absolute flow. Hence, 


Tplosa 


= \p E — 'pm = 


2 (]II 1033 

U? 



( 10 ) 


where H\ 0%a = H E — H M is the head loss due to friction in relative flow. 

For any given inducer, ip E — \p m can be determined from the measured 
stagnation pressures and tangential velocities at the exit, and, knowing 
the inducer geometry, values of X can be derived at any radii. The loss 
factor derived from various measurements carried out at NASA Lewis 
Research Center, M.I.T. Gas Turbine Laboratory, and TRW Cleveland 
are plotted in figure 8. The sources from which these data were collected 
and details of the inducer configurations and flow parameters used are 
given in table II. The radial distribution of Blasius loss coefficient showed 
a systematic trend from midradius to tip ; i.e., an increase in frictional 
losses with decreasing flow coefficient. Dependency of energy losses on 
hub/tip ratio is also evident from figure 8. Larger secondary losses asso- 
ciated with low hub-tip ratios are well known. A new friction factor X« 
is defined and all the data are replotted to derive the value of this new 
friction factor. 
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Figure 8. — Radial variation of 
Blausius friction coefficient (X) 
evaluated from equation (10) 
for various inducer configura- 
tions ( table II). 




Table II . — Details of the Various Inducer Configurations and Flow Parameters Used in Deriving the Friction Factors X and ( Figures 8 and 9) 
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1 All inducers listed were tested in water. 
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2f/// loss ^ Rht 1 L /II V 


( 11 ) 


whore Rht is the hub/tip ratio and H\ oas = APn/pg. 

In these calculations, local values of blade length, hydraulic diameter, 
blade mean angle (arc tan (tan ft + tanft)/2), design flow coefficient, 
and measured local velocity (IT) at the outlet are used to derive the new 
friction factor (A ff ) for the various inducer channels (fig. 9). It is not 
certain whether such a linear dependency on fir, /IT* exists, but the data 
plotted in figure 9 seem to indicate that a relationship such as equation 
(11) is closely obeyed by all the inducer data plotted. 

It should be remarked here that these loss correlations include mixing 
losses downstream of the; trailing edge, since the measurements on which 
these loss correlations are based are taken for downstream of the inducer 
and not at the trailing edge. 

It is suggested that a relationship such as the one discussed above be 
used with values of distribution shown in figure 9 for the flow coefficient 
c/> ranging from 0.065 to 0.2. The largest scatter in the data is due to the 
M.I.T. inducer which was tested at 0 = 0.2. The secondary losses at the 
annulus and hub walls tend to be as important as secondary flow inside 
the blade boundary layer in such a case. 



Figure 9. — Radial variation of 
the modified friction loss co- 
efficients (X R ) calculated from 
equation (11) for various in- 
ducer configurations ( table II). 
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Circumferential Average of Radial Equilibrium Equation 

The circumferential averaging of various terms in equation (7) is 
carried out under the following assumptions (fig. 7). 

(1) The passages have fully developed turbulent flow and the flow- 
properties are symmetrical about the midpassage. 

(2) The velocity component parallel to the blade passage is given by 

W x — W m v ,n (12) 

where ll’ m is the velocity component parallel to the blade at midpassage, 
i) = y/S, where y is the distance measured normal to the blade, 2>S is the 
width of the channel (fig. 7) and x is the coordinate parallel to the blade 
and lying on a cylindrical surface. 

(3) The radial velocity component W r is given by reference 12, 

W r = W T tan a ( 1 - r,) 2 = W m tan a ( v) 1/7 ( 1 - r,) 2 (13) 

where a is the angle between the direction of the resultant velocity near 
the wall and W m or angle between the limiting streamline and x direction 
(fig. 7) . Note that W m and a are functions of r and x only. The average of 
any quantity T will bo denoted by a bar and is defined by 

T = \fTdy (14) 

*5 J o 

where 2S is the width of the blade passage. 


Circumferential Average of Term 1 in Equation (7) 

Since W x is the velocity along x, where x is an intermediate variable 
and r, 6 are the independent variables, 


dW r TI7 dW 

x = W z 

dx dz 


~+w, 


fllTr 

rdd 


(15) 


Now, from equation (13), 


dW r n( ,j 
dx 


dW m Trr , da 

— tan a+ H m sec- a — 
dx dx 


( 16 ) 


Hence, the circumferential average of Term 1 in equation (7) is 


We dWr dWr dWr 
r dd + Z dz x dx 


= Wrr 


tan a+W m sec 2 — \ f T n (l — y) 2 dy 
dx dx/ Jr, 
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= 0.210 W m 



tan a-\-W m sec 2 a 


3a\ 

dx) 


(17) 


Circumferential Average of Terms 2, 3, and 4 in Equation (7) 


Since it is assumed that the fluid is smoothly guided through the 
passage, the circumferential average of Term 2 in equation (7) is 


W 


aw, 

dr 


= f (W e ) m (v) in (v) ind -^~d v =1.015We d ™ 

J ft 


dr 


dr 


(18) 


since We = \(We)m for one-seventh profile. 

Similarly, 

WV= f We 2 dr, = (We) J [ r n d v = 1.015W? 

•'o ■'o 

and 

2nWe = 2QWe 


(19) 

(20) 


Circumferential Average and Radial Variation of the Frictional Losses 
in Equation (7) 


Since the friction losses derived by the author (eq. 11) represent the 
average losses across the passage, the circumferential average of the fifth 
term in equation (7) is 


d/APo\d/ Rht 1 L 
dr\ p ) dr V* <t> R N m d h 2 ) 


( 21 ) 


Hence, the circumferential average of the radial equilibrium equation now 
reads, after substituting equations (17), (18), (19), (20), and (21) in 
equation (7), 


— 0.210 W m ( tan . 


dW„ 


dx 


3 \ 

-\-W m sec 2 a. — ) + 1.015TFe 


dx/ 


yh+i.ms’t*’ 

dr r 


- 21211 e~\~~ 


R h t L 1T 2 \_ 0 


dr \R N Ui <f> dh 2 ) 


( 22 ) 


Since We = J (We) m = : gW m sin /3, the final circumferentially averaged radial 
equilibrium equation valid for inducers operating at low flow coefficient is 

-0.27 5 
sin 2 /3 


W, 


( 


dWe da\ 311 e 

tana bJF«sec 2 a — 1+1.01511^ 

dx dx) dr 


We 2 

+ 1.015 — — 2SlWe 
r 


A+^ L .!^) = o (23) 
dr\R N l/i 4> d h 2 1 
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The solution of this circumferentially averaged radial equilibrium 
equation to derive We or Ve distribution along the radius is still intractable 
unless suitable assumptions are made for dWe/dx and a. The experiment 
is the only source from which the values of a can be derived. The limiting 
streamline angle a is predictable in the case of a single rotating helical 
blade (ref. 5), but the values of a measured on inducer blade surfaces are 
found to be higher than those on a single helical blade (refs. 21 and 24). 
The evaluation of dWe/dx could be based on either a linear variation of 
We along x or the inviscid solution. Once these quantities arc evaluated, 
We or Ve can be predicted from equation (23) using the author’s correla- 
tion for friction loss coefficient X« (fig. 9). 

The method of predicting We or Ve for any given inducer geometry and 
flow coefficient is described in the next section. 


EXPERIMENTAL RESULTS AND DISCUSSION 

The four-bladed inducer described in the section on the design of the 
test inducer is used for the experimental investigations of the flow char- 
acteristics. The test medium is air and all the measurements are carried 
out at a Reynolds number (based on tip radius) of 6.6 X10 5 , which 
corresponds to 450 rpm. A throttle valve, which can be seen in figure 2, is 
adjusted to obtain a flow coefficient <t> of 0.065. 

The measuring stations are shown in figure 10. The exit flows were 
measured at stations 3, 4, 5, 6, and 7, using a stagnation probe and a tuft 
probe. Hot wire measurements were taken very close to the trailing edge 
(station 3A) to derive the average values of radial, axial, and tangential 
velocity components. 

The outlet air angle's (of the absolute flow) were measured by means of 
a very thin tuft probe and a sighting scope of the type used by Smith 
(ref. 22) . Fluctuations in absolute air angles were found to be consider- 
able, especially near the hub locations. The algebraic mean of the maxi- 
mum and minimum air angles, measured by the tuft probe, was used for 
calculating the axial and tangential velocity components. In spite of slight 
inaccuracies which are discussed by McCafferty (ref. 20) , much informa- 
tion can be gained from these stationary probe measurements. Since the 
flow at the exit was highly three-dimensional, the static pressures were 
measured at the hub and annulus walls and the radial distribution of 
static pressure coefficient was derived by computing the pressure gradient 
near the wall locations from the simplified radial equilibrium equation. 

Since an accurate knowledge of the nature and magnitude of radial 
velocities is essential for the development of an adequate theory for pre- 
dicting the flow, an attempt was made to derive the average values of all 
three velocity components at the exit (station 3A) using a combination 
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Figure 10. — Location of the flow-measuring stations. 


of an A'-configuration hot-wire probe plus a single sensor probe normal to 
the A'-configuration. The probes were aligned in the three coordinate 
directions such that they sense the following velocity components. 


V e >= VVe 2 +V* 

(24) 

V er = VV7+ Vr 2 

(25) 

v rz = Vv7+V7Tkw? 

(26) 


Since the air angles of the absolute flow are very large, a correction to 
the cosine law was used for the wire in the 6 direction. This appears in the 
form of a correction factor as shown in equation (26). The value of K was 
derived from Schwarz and Friehe’s correlations (ref. 23). For the length/ 
diameter of the wire used in these experiments, the value of K was found 
to be 0.26. 

The oscilloscope traces of the signal were photographed 4 to derive the 
blade-to-blade variation of the voltages. The hot-wire calibration curves 
and equations (24), (25), and (26) were then used to derive the blade- 
to-blade variation of V r , V e , and V z . 


4 The method of data processing has since been improved as indicated in reference 4. 
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Tangential Velocities at Exit 

The absolute tangential velocities derived from pressure and angle 
measurements are shown plotted in figure 11 for various axial stations. 
The measured tangential velocities are substantially higher than the 
design values except near the hub. These results are consistent with the 
measurements reported by Sandercock et al. (ref. 17) and Soltis et al. 
(ref. 18). The real fluid effects combined with large radial flows that exist 
inside the blade passages are responsible for such a large increase in 
tangential velocity or head rise. Furthermore, it is also evident from 
figure 11 that there is a substantial change in tangential momentum as 
the flow proceeds downstream. Far downstream of the inducer, a “forced 
vortex” type of distribution exists. Such changes are due to mass flow 
redistribution that takes place downstream of the inducer and the mixing 
losses due to blade wake diffusion. The extent of mass flow redistribution 
can be seen from figure 1G. These measurements indicate the desirability 
of having the inducer as an integral part of the main pump and also 
point out the need to take measurements close to the trailing edge in order 
to accurately assess the extent of real fluid effects inside the blade passages. 
It is quite clear that the inducer behaves like a drag or a friction pump, 
especially near the tip, where the head rise is imparted to the fluid by 
purely frictional effects. 


Figure 11 . — Radial variation of 
absolute tangential velocity at 
various axial locations. 
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Even though an accurate prediction of the tangential velocity dis- 
tribution is possible only by a thorough understanding of the three- 
dimensional boundary layer characteristics inside the passage, a qualita- 
tive estimate of the frictional effects can be made from an analysis of the 
type indicated in the preceding section (Approximate Theoretical 
Analysis) . 

The form of equation used for predicting the tangential velocities at 
the exit of the inducer is obtained by further simplification of equation 
(23) (by expressing We as the only dependent variable) under the 
following assumptions: 

(1) lt'« varies linearly with x from leading to trailing edge. 

(2) Angle a between the limiting streamline and W m varies linearly 
from leading edge to trailing edge, being zero at the leading edge. This 
trend is confirmed from the limiting streamline angle measurements 
reported in reference 21. Hence, 


dWe We — fir We — Or 
dx L 4.5 r 


(27) 


where 

We = relative tangential velocity at exit 

fir = relative tangential velocity at inlet 

L = chord length of the blade (=4.5r for Penn State inducer) 

and 


da a — 0 

dx 4.5 r 


(28) 


where a is the angle of the limiting streamline near the trailing edge. Use 
of equations (27) and (28) would reduce the first term in equation (23) 
to 


where 


-We 2 

r 


f(r)+2VT,fi$(r) 


f(r) = 


0.061 
sin 2 /3 


(tan a+a sec 2 a) 


(29) 


(30) 


Hr) 


0.03 tan a 
sin 2 /3 


(31) 


The friction loss term in equation (23) can be further simplified for the 
Penn State inducer to give 
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d_ 

dr 




4> J>V /4 4 2 / dr 


0 1 1 ft 

= F'(r)Wf+2F(r)We (32) 

dr 


where 


F (r) = 


1 

2 


/V 74 


Rkt L 1 

<f> d h sin 2 ft 


(33) 


and \ R (fig. 9), f h , (oq. (9)), length of the blade (L), and sin 2 fi 

(fig. 6) are all functions of radius at any given axial station. The variation 
of F{r) with radius at the location of the trailing edge is plotted in figure 
12. Hence, the final form of radial equilibrium equation used for pre- 
dicting the relative tangential velocity at the exit of the Penn State 
inducer is obtained by substituting equations (29) and (32) in equation 
(23). 


Figure 12. — Radial variation of 
function F(r) ( eq . (33)) for 
Penn State inducer. 
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F,— (1.015+2F(r))+ — (1.015+rF , (r)-r(r))=2aF,(l-«(r)) 
dr r 

(34) 

where the functions F(r), f(r), and £(r) are defined by equations (33), 
(30), and (31), and F'(r) is the derivative of F(r) plotted in figure 12. 

Equation (34) is used to predict the radial variation of mean relative 
tangential velocity (We) at the trailing edge. In evaluating the functions 
f(r) and £(r), information obtained by flow visualization experiments 
reported in reference 24 is used. Near the hub, the value of a was found 
to be nearly 60° on the suction surface and 30° on the pressure surface of 
the blade. Basing our assumption on this, a is assumed to vary linearly 
from 45° at the hub to zero at the tip (fig. 13) . This enabled evaluation of 
the functions f(r) and £(r) (eqs. (30) and (31)) in equation (34). The 
variations of these functions with radii are plotted in figure 13. 

The predicted values of relative tangential velocity at the trailing edge 
obtained from equation (34), using estimated values of F(r), F'(r), f (r), 
and £(r) (figs. 12 and 13), are plotted in figure 14. The boundary condi- 
tion assumed for the solution of the differential equation (34) is that the 
tangential velocity (We) at the midradius is the same as that of design. 
The predicted values are in good agreement with the values derived from 
stationary probe measurements at station 5. 5 

The effect of neglecting the inertial terms due to radial velocity gradients 
(i.e., functions f(r) and £(r) in equation (34)) is also shown in figure 14. 
The inertial terms due to radial velocity gradients in the radial equilib- 
rium equation seem to have very little effect on the predicted values of 
We- Even though the acceleration components due to radial velocity 
gradients are small, their influence in increasing the frictional losses 
cannot be ignored. Hence, in conclusion, it can be said that the accelera- 
tion component due to radial velocity gradients can easily be neglected in 
predicting the relative tangential velocity distribution from the radial 
equilibrium equation. But the influence of radial velocities on axial 
velocity profile at the exit cannot be ignored. 

In deriving equation (34) from the generalized radial equilibrium 
equation (eq. (4)), it was assumed that the term W z (dW z /dr) is very 
small. Its influence on the predicted values of Wo is evaluated by using the 


5 It should be remarked here that the loss coefficient (\h) derived in this paper is 
based on the measurements carried out at one-seventh to one-ninth blade chord (refs. 
15, 17, and 18) downstream of the trailing edge. Hence, these loss correlations include 
mixing losses. To be consistent, only the measurements carried out at station 5, which 
approximately corresponds to the measuring stations of references 15, 17, and 18, is 
compared with the predicted values. 
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We 

Ut 


Figure 14. — Predicted and meas- 
ured relative tangential velocity. 
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experimental values of W z (dW z /dr) (fig. 16) and the following modified 
equation 

a it a tt+ dw z 

W,— (1.015+2 F(r))+— (1.015+rf , '(r)-f(r))+H'. — 
dr r dr 

= 2flFi(l-f(r)) (35) 

The effect of neglecting the axial velocity gradient term in equation 
(34) is evident from figure 14. The effect of this term on the predicted 
values of We seems to be appreciable only near the hub (where the axial 
velocity gradients are large) and negligibly small at other radial positions. 

The predicted tangential velocity components of the absolute flow 
(U—We) are plotted and compared with experimental values, derived 
from stationary probe measurements 6 at station 5, in figure 15. The agree- 
ment between the predicted and experimental values is good. The inclu- 
sion of the axial velocity gradient term W z {dW z /dr) in equation (34) 
seems to bring the predicted values of Ve closer to the experimental 
values only near the hub. 



Figure 15. — Predicted and meas- 
ured values of absolute tangen- 
tial velocity. 


9 See footnote 5. 
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The effect of neglecting the radial velocity gradients f (r) and £(r) in 
equation (34) is shown in figure 15. It is evident that the friction loss 
term in the radial equilibrium equation (cq. (7)) has a very appreciable 
effect in changing the tangential velocity gradient in the radial direction. 

The tangential velocity distribution derived from hot-wire probe and 
stationary probe measurements, carried out very close to the trailing 
edge (station 3A), is compared in figure 17. There is a slight discrepancy 
between the two measurements. This may be due to error in measuring the 
fluctuating (with reference to stationary probe) stagnation pressure close 
to the trailing edge. 

Axial Velocities at the Exit 

The measured radial distribution of axial velocities at various axial 
stations is plotted in figure 1G. 

Immediately downstream of the inducer (stations 3 and 4), the mini- 
mum velocity occurs at the midradius. This seems to suggest that radial 
migration inside the blade boundary layer is very large from midradius to 
tip as compared to that from hub to midradius. At station 3, where the 
hub flow is still under the influence of blade, the blockage effects are large. 
Sandercock ct al. (ref. 17) and Soltis et al. (ref. 18) concluded earlier that 
the inducers operating at low flow coefficients have separated flow near 
the hub. These conclusions are based on measurements taken far down- 
stream. The present investigation seems to suggest that such flow separa- 
tion takes place downstream of the inducer due to considerable flow re- 
distribution. It should be remarked here that large radial velocities exist 


Figure 16. — Radial variation of 
axial velocity at various axial 
locations. 
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in the blade wakes of the inducer. This accounts for large changes in axial 
and tangential velocities and stagnation pressures measured downstream 
(figs. 11, 1G, and 18). 

Far downstream of the inducer (stations 5, G, and 7), the minimum 
axial velocity occurs near the hub, and in fact shows a tendency to sepa- 
rate at station 5. Downstream of the inducer, the radial velocities inside 
th(' blade wakes are outward near the hub, whereas the wake interactions 
near the tip produce an inward flow. The net effect is a migration of the 
mass flow toward the 1 midradius. This effect can be clearly soon at stations 
5 and 6. Thus, there is a substantial improvement in axial velocity distri- 
bution near the midradius and a deterioration near the hub. Furthermore, 
the extent of back flow is observed to grow continuously as the flow pro- 
ceeds downstream and occupies nearly 20 percent of the annulus far down- 
stream (station 7, fig. 1G). In a pumping unit where the inducer is an inte- 
gral part of the main pump, the back flow may be absent if aerodynam- 
ically designed blade sections are used. It is not clear whether this is true 
for flat plate inducers, since measurements close to the trailing edge are 
not available for such units. Inside 1 the blade passages, the perturbation in 
axial velocity is due to : 

( 1 ) Radially varying blade thickness and hence blockage 

(2) Radial flow inside the blade boundary layers. This is by far the 
most important cause of the change in axial velocities from hub to tip. 

The perturbations in downstream flow are predominantly due to radial 
velocities in the three-dimensional wakes. 

The prediction of axial velocity distribution should take into con- 
sideration the detail flow structure, including the radial velocities inside 
blade boundary layers, inside such inducers. The boundary layer in a 
rotating helical channel, currently under investigation at Pennsylvania 
State University, should provide the urgently needed information for the 
accurate prediction of axial velocity distribution. Thus, considerably 
more theoretical and experimental work is necessary before the observed 
flow can be regarded as rationally explained. 

It should be emphasized here that the flow redistribution downstream 
of the inducer (caused by the radial velocities in the three-dimensional 
blade wakes) is considerable and extreme caution should be exercised in 
evaluating the performance of an inducer from the measurements obtained 
far downstream. 

The average values of the axial velocities derived from the hot-wire 
probe measurements very close to the trailing edge are shown plotted and 
compared with stationary probe measurements in figure 17. This seems to 
confirm the trend in axial velocity discussed earlier; i.e., very low axial 
velocities at midradius. 
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Figure 17. — Comparison of the absolute velocity components measured by the hot wire and 
stationary probe at station 3A. 

Radial Velocities at Exit 

The radial distribution of average radial velocities (| V r |) derived from 
the hot-wire measurements is plotted in figure 17. 7 

It should be remarked here that the hot-wire probe senses only the 
magnitude and not the direction of the radial velocity. Near the inducer 
tip, radial inward flow exists due to boundary layer interactions. Thus, 
the algebraic average of radial velocities across the passage is likely to be 
different from | V r \, especially near the tip. 

It is evident from figure 17 that the radial velocities are of the same 
order of magnitude as axial velocities, with values increasing monotoni- 
cally toward the tip. The maximum average radial velocity occurs a small 
distance away from the tip (fig. 17) and this location corresponds to the 
maximum inward flow region. 

These measurements confirm the earlier conclusion that radial velocities 
in such low flow coefficient and high solidity inducers are large. For the 
eventual flow analysis, it is essential to know the detail blade-to-blade 
variation of radial velocities including the extent of radial inward flow. 


7 Discrepancy between the hot-wire results reported in reference 21 (unpublished) 
and this paper is due to the fact that earlier results are based on K = 0 in equation (26). 
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Stagnation Pressure Coefficient at Exit 

The radial variation of stagnation head coefficient at various axial 
stations is shown plotted in figure 18 and compared with design values. 
These plots are typical of inducers tested earlier (refs. 17 and 18). The 
pressure coefficients near the tip are two to three times the corresponding 
design values. Strong real fluid effects from midradius to tip arc thus 
evident. The performance of the inducer is similar to a “shear force pump” 
or “drag pump,” where shear forces are utilized to exchange kinetic 
energy between fluid and solid bodies. Thus, the sections near the tip 
operate predominantly as a shear force pump, whereas the pressure rise 
near the hub is largely due to flow turning effects. 

The stagnation head decreases continuously as the flow proceeds down- 
stream of the trailing edge. Major changes occur between stations 3 and 
4, especially near the hub. Comparing the corresponding axial velocities 
at stations 3 and 4 in figure 16, it is evident that the major effects come 
from the axial velocity changes. The wake mixing losses also contribute 
to the change in stagnation head. These effects are dominant from mid- 
radius to tip. The tip sections experience considerable pressure losses as 
the flow proceeds downstream despite the fact that axial velocity changes 
(fig. 16) are small. The presence of large mixing losses is thus evident in 
these regions. 

An attempt is made to estimate the stagnation pressure or head rise 
coefficient and losses on the basis of the analysis indicated in the preceding 
section (Approximate Theoretical Analysis). The knowledge of the 



Fiuuue 18. — Radial variation of 
absolute stagnation pressure co- 
efficient ( i /' t ) at various axial 
locations. 
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relative velocities derived from the theory outlined earlier (see Tangential 
Velocities at Exit and fig. 15) enables an estimate of the energy losses to 
be made from equation (11). The loss coefficients so derived are shown 
compared with the measured losses in figure 19 for station 5. 8 Reasonably 
good agreement between the estimated and measured losses confirms the 
validity of an empirical correlation such as equation (11) for estimating 
the energy losses associated with friction. 

Since the theoretical values of Euler head {\p E = 2RVe/U t ) are known 
from the predicted values of Ve (fig. 14), the actual head rise can be esti- 
mated from the equation \J/ t = 'Pe — 'I' i OS8 . 

The value of \p T thus estimated is shown plotted and compared with 
measured values at station 5 in figure 19. Here again, the agreement is 
reasonably good and large gradients in stagnation head near the tip arc 
estimated qualitatively. ,, ' 

Static Pressure Coefficient at the Exit 

The static pressure coefficients derived at various axial locations 
arc plotted and compared with design values in figure 20. 



Figure 19 . — Estimated and measured friction losses and stagnation pressure coefficient 

0/t). 


8 See footnote 5. 
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Figure 20. — Radial variation of 
V's at various axial locations. 


Near the tip, the measured pressure gradients are found to be large 
compared to those of the design. This is in conformity with the radial 
equilibrium requirements since tangential velocities at these locations are 
very large. It .should be remarked here that the hub flow at station 3 is 
still under the influence of the rotor and hence the increase in \ps observed 
between stations 3 and 4. There is appreciable loss in static pressures as 
the flow proceeds downstream. 

SUMMARY AND CONCLUSIONS 

(1) A new friction loss coefficient applicable to inducers operating in 
the range of flow coefficients <f> = 0.0(55 to <t> = 0.2 is defined and derived 
from the inducer data available in the literature. This empirical friction 
loss coefficient is found to increase exponentially toward the blade tip. 
For the Penn State inducer, the radial variation of frictional losses esti- 
mated from this newly derived empirical loss coefficient agrees closely with 
the measured values. 

(2) A circumferentially averaged radial equilibrium equation is used 
to predict the relative and absolute tangential velocities. The analysis is 
based on suitable assumptions for the radial and mainflow velocity profiles 
(based on the existing three-dimensional turbulent boundary layer data 
available) and the loss coefficients explained above. The agreement 
between theory and experiment is reasonably good. Hence it is evident 
that, if the frictional effects are known either empirically or analytically, 
the flow properties at the exit of the inducer can be predicted quite 
accurately. The stagnation head rise estimated from the predicted 
tangential velocities and the stagnation pressure losses (item (1)) agree 
qualitatively with the measured values. 
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(3) The test inducer, designed approximately for uniform head 
distribution over its discharge area (assuming ideal flow), actually pro- 
duces a nonuniform head. Near the tip, the actual head of the absolute 
flow was found to be two to three times that at hub and midradius. This 
nonuniform head distribution can be explained qualitatively by real fluid 
effects and is in agreement with the observations of other investigators. 
The stagnation pressures are found to decrease continuously as the flow 
proceeds downstream. The wake mixing losses combined with large change 
in axial velocity distribution are the major causes of this. 

(4) The expected radial motions within the blade passages have been 
qualitatively confirmed by hot-wire measurements and appear to be quite 
strong at all radii. The radial velocities are found to be of the same order 
of magnitude as axial velocities and increase continuously towards the tip. 

(5) Measurements carried out several stations downstream indicate 
that the axial velocity profiles undergo marked changes as the flow pro- 
ceeds downstream. The back flow region develops after the flow has left 
the rotating hub and grows continuously as the flow proceeds downstream. 
The minimum velocity occurs near the; midradius, immediately down- 
stream of the trailing edge, whereas, far downstream, the hub locations 
have minimum axial velocity. 

(6) The absolute tangential velocities are substantially higher than 
the design values, except near the hub. These are consistent with the 
stagnation pressure* rise observed. There is a large change; in tangential 
velocities as the flow proceeds downstream. This is caused by mass flow 
redistribution and mixing losses that occur downstre>am of the inducer. 

(7) The measures! static pressure; gradients are found to be large 
compared to those of design, especially near the tip. There is appreciable 
loss in static pressures in the tip flow as the flow proceeds downstream. 
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LIST OF SYMBOLS 

dh Hydraulic mean diameter (equation (9)) 

F Frictional force per unit mass 

F{r), t(r), £( r ) Functions defined in equations (33), (30), and (31) 
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Stagnation head rise (in absolute flow) 

Head rise derived from Euler’s Equation ( UVe/g ) 
H E -H m 

Measured stagnation head rise (in absolute flow) 

Net positive suction head, ft 

Static head 

Blade chord length 

Rotative speed, rpm 

Number of blades 

Static pressure 

Stagnation pressure 

Stagnation pressure loss due to friction 

Flow rate, gpm 

r/r t 

Hub/tip ratio 

Reynolds number, W dh/v 

Radial, tangential, and axial coordinates 

Width of the passage measured normal to blade surface 

Suction specific speed, N y/ Q/H sv 3li 

Blade spacing 

Blade speed 

Resultant absolute velocity 
Resultant relative velocity 

Coordinates parallel to and perpendicular to blade 
camber line and lying on a cylindrical surface (fig. 7) 
Angle between limiting streamline and x direction 

(fig- 7) 

Blade angle measured from axial direction 
Air angle measured from axial direction 
y/S 

Distance between the mean streamline and camber line 
normal to blade chord 

Blasius friction coefficient defined in equation (8) 
Friction coefficient for a rotating channel defined by 
equation (11) 

Kinematic viscosity 
Density 

Flow coefficient, lf\/flr, 

2 gH E /U? 

2gH i OS3 / U i 2 

2gh/U t 2 (static head rise coefficient) 

Angular velocity 

Stagnation head coefficient, 2 gH/ U ? 
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Subscripts 


h 

m 

RMS 
r, d,z 
t 


x, y 

1 

2 


Hub 

Values at midpassage 
Root-mean-square value 
Components along r, 6, z directions 
Tip 

Components along x, y directions (fig. 7) 

Inlet 

Outlet 


Superscript 

Averaged over the passage 


REFERENCES 


1. Acosta, A. J., An Experimental Study of Cavitating Inducers. Proc. Second 

Symposium on Naval H ydrodynamics (Washington, D.C.), August 25, 1958. 

2. Arcand, L., The Performance of Two Axial Flow Water Jet Pumps. ASME Sym- 

posium on Pumping Machinery for Marine Propulsion, May 1968. 

3. Poncet, A., and B. Lakshminarayana, Investigations and Analysis of Flow 

Phenomona in A xml Flow Inducers. NASA CR 107267, October 1969. 

4. Poncet, A., and B. Lakshminarayana, Investigations of Three-Dimensional Flow 

Characteristics in a Three Bladed Rocket Pump Inducer. NASA CR-2290, 1973. 

5. Lakshminarayana, B., A. Jabbari, and H. Yamaoka, Turbulent Boundary 

Layer on a Rotating Helical Blade. ,/. Fluid Mech., Vol. 51, p. 547, 1972. 

6. Wisucenus, G. F., Fluid Mechanics of Turbomachinery. Dover, Inc., 1965, 

Chapter 29. 

7. Wislicenus, G. F., and B. Lakshminarayana, Design of a Test Inducer. Progress 

Report, NASA Contract NSG 537, March 1964. 

8. Ackeret, Contribution to the Layout of Closely Spaced Cascade of Vanes. Schwei- 

zer Bauzeitung, Vol. 120, No. 9, August 1942. 

9. Cooper, P., and H. Bosch, Three-Dimensional Analysis of Inducer Fluid Flow. 

NASA CR 54836, TRW ER 6673A, February 1966. 

10. Rtjden, P., Investigations of Single-Stage Fans. NACA TM 1062, April 1944. 

11. Smith, L. H., The Radial Equilibrium Equation of Turbomachinery. Trans. 

ASME, Ser. H, Vol. 88, 1966, pp. 1-12. 

12. Mager, A., Generalization of Boundary Layer Momentum Integral Equation to 

Three-Dimensional Flows Including Those of Rotating System. NACA TR 1067, 
1952. 

13. Spannhake, W., Proc. Nat. Conf. Irul. Hydraulics, Vol. 8, October 1944. 

14. Ludwieg, H., Ing. Arch., Vol. 19, No. 296, 1951. 

15. Mullan, P. J., An Investigation of Cavitating Inducers for Turbopumps. M.I.T., 

Gas Turbine Lab Report 53, May 1959. 

16. Bosch, H. B., et al., Advanced Inducer Study. ER-5288, NASA Document N63- 

21124, TRW, Inc., May 1963. 

17. Sandercock, D. M., et ad., Cavitation and Noncavitation Performance of an 80.6° 

Flat Plate Inducer. NACA TN D1439, November 1962. 



INVESTIGATION OF FLOW DISTRIBUTION IN ROCKET PU.MP INDUCERS 723 


18. Soltis, R. F., et al., Investigation of the Performance of a 78° Flat Plate Helical 

Inducer. NACA TN D1170, March 1962. 

19. Vavra, M. H., Aerolhermodynamics and Fluid Flow in Turbomachines. John 

Wiley and Sons, Inc. (New York), 1960, p. 123. 

20. McCafferty, H. G., Errors in Measuring the Fluctuating Flow at the Discharge of 

an Inducer. M.S. thesis, Dept. Aerospace Eng., Penn State U., June 1967. 

21. Lakshminarayana, B., Investigations and Analysis of Flow in Axial Flow Inducer. 

NASA CR 103291, June 1969. 

22. Smith, L. H., Three-Dimensional Floiv in Axial Flow Turbomachinery. Johns 

Hopkins U., Report 1-14, November 1953. 

23. Schwarz, W. H., and C. A. Friehe, Deviations from Cosine Law for Yawed Cy- 

lindrical Anemometer Sensors. ASME Paper C8-WA/APM-16. 

24. Lakshminarayana, B., Visualization Study of Flow in Axial Flow Inducer. ./. 

Basic Eng., Dec. 1972, p. 777-787. 


724 PUMPING MACHINERY FOR AEROSPACE APPLICATIONS 


DISCUSSION 


G. F. WISLICENUS (Tucson, Arizona) : I am very happy to discuss 
this paper, as I was connected for a considerable length of time with the 
work on which this paper reports. 

Our friends from NASA will recall that we had no illusions about the 
inherent difficulties of this investigation at the time it was proposed. As I 
recall it, we stated in this proposal that real-flow effects were expected to 
extend over the entire flow field. In other words, we expected to find a 
mess, and this expectation was certainly fulfilled. The objective we had in 
mind was not only to learn something about rocket pump inducers, but 
equally to learn how to deal with such a complex and theoretically almost 
hopeless flow problem. It may be understandable why an old engineer like 
myself might want to get involved in such a problem. But it commands 
my respect when a younger engineer and scientist like Dr. Lakshminara- 
yana is willing to struggle for years with this type of investigation, and 
even more so if he finally comes up with results that seem to make sense. 
I can assure you that while I was still at Penn State I found hardly any 
time to contribute to this often frustrating work. And now I am com- 
pletely out of it, which has given Dr. Lakshminarayana the opportunity 
to pursue his investigation in a more undisturbed fashion. 

I do not feel qualified to discuss the author’s theoretical investigations 
except for the observation that some of the comparisons between his 
theoretical and experimental results impress me as much closer than I 
personally would have expected. But if it is really true, as stated in con- 
clusion (2), that frictional effects in inducers are now predictable, then it 
is our obligation to put this knowledge in a form that can be used by the 
practical designer. After all, if friction effects double or triple the head 
generated in the tip regions, then their consideration is a plain necessity 
in the design process. 

There is only one point I should like to discuss in a little more detail; 
i.e., that according to figures 11 and 18 the angular momentum of the 
flow seems to decrease with increasing distance from the impeller. The 
statement that this is due to “mass flow redistribution” might not be a 
sufficient explanation for some old-fashioned engineers like myself. It 
seems to me that this situation could be cleared up by plotting the cir- 
cumferential velocity component (To) against the stream function of the 
(circumferentially averaged) meridional flow, as the local angular 
momentum transport is, after all, proportional to the local mass flow. If 
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such a plot still shows a decrease in angular momentum, this decrease 
must be explained by a circumferentially nonuniform mass and momen- 
tum exchange between the positive flow region and the backflow region. 

In passing, I should like to ask whether the axial velocity curves shown 
in figure 16 satisfy the condition of continuity. To check this, it might be 
interesting to plot the same curves against r 2 , considering the backflow 
areas as negative. It illustrates the measuring problem before us if we 
mention that under such complex conditions of flow measurement any 
check on the validity on the results obtained is most welcome. 

It worries me just a little that the simple throttle valve which I sug- 
gested to use might have an upstream effect that could influence the 
separation at the hub downstream of the inducer. This possibility could 
be checked with relative ease. 

J. E. CROUSE (NASA Lewis Research Center) : The author and 
investigators associated with this project should be congratulated on the 
depth of their research work. It is encouraging to see some “in depth” 
studies among the large number of papers published these days. 

In the paper, considerable effort was made to construct models of 
velocity profiles and loss mechanisms. The data essentially confirms that 
suspected significant radial profiles near the blade surfaces do indeed 
exist. One of the more interesting results of the work is the small effect 
these radial components have on the predicted energy addition as evi- 
denced by the V>/U t profiles of figure 15 with and without the radial 
velocity terms. 

The important empirical term in the author’s flow prediction equation 
is the term representing the loss gradient. Especially for low flow co- 
efficient turbomachines, such as inducers, a good flow prediction method 
is highly dependent on an accurate representation of the loss gradient. 
The source of the losses is not important for the prediction method itself, 
but knowledge of the source of the losses can be rather important in 
setting up a reasonable loss model. The author chose to correlate all the 
losses associated with the flow through an inducer into a pipe-type friction 
parameter for the relative flow channel. 

The experimental data of figure 9 for several rotors, according to the 
author’s suggested correlation, is indeed somewhat better than that of 
figure 8, which is the same data by another correlation. The better correla- 
tion of loss data should be considered one of the most important results 
of this work. 

J. H. HORLOCK (Cambridge University) : To one who has attempted 
to predict the averaged flow in conventional axial flow turbomachines, 
this is a fascinating paper in that the method of attack has to be entirely 
different for several reasons. 
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(1) The axial velocity term in the radial equation of motion 
W z (dW z /dr) is neglected because of the high stagger of the blade. 

(2) Prediction of the tangential velocity distribution TT«(r) instead 
of W z (r) is attempted. 

(3) Variations of velocities across the pitch due to viscous effects 
are important. 

(4) The radial loss distribution is dominant (although this is also 
important in axial-flow compressors at low flow coefficients). 

I have followed B. Lakshminarayana’s work with interest and discussed 
it with him. But there are still several aspects of the investigation that are 
not clear to me. 

(1) On what basis is it justifiable to neglect dW z /dr compared with 

dWr/dz? 

(2) In the region near the blade in the outer half of the rotor, I think 
the magnitude of the radial body force due to viscous loss should be 
estimated. It is of order of magnitude [tan a d(AP<i/p)/d.r] and may 
not be negligible. 

(3) Why is the loss correlation improved by including the flow coeffi- 
cient </>? I would expect it to be related to (4> — <t> design) but not simply as 
inverse of <f>. 

(4) The change in Ve(r) with r (fig. 11) is mysterious. Is it possible 
to measure the torque on the rotor and compare it with the change of 
tangential velocity across the rotor, in order to see if there are errors in 
the measurement of the latter. 

Undoubtedly this work will be of use to designers of other turbo- 
machines, especially in the prediction of performance at low <f>. Has the 
author any recommendations on how the design of the inducer may be 
modified to allow for the observations of this paper? Can the large in- 
crease in loss with radius be avoided, and is the operation of the outer 
part as a drag pump inevitable? Or should the inner part be designed 
inviscidly and the outer part as a drag pump? 

P. COOPER (Case Western University and TRW Inc.): Dr. 
Lakshminarayana’s work makes one aware of the strong three-dimensional 
effects that serve to alter the expected performance of an inducer. He has 
made us aware of heretofore unexplained radial flows and other effects of 
the highly frictional flow that exists in an inducer. That he has been able 
to predict the outlet tangential velocity distribution by using his new loss 
scheme is especially encouraging. 

I would like to inquire about the flow process connected with his pre- 
dicted shift in this tangential velocity distribution from that of design. 
The author suggests that the flow action near the inducer tip is like that 
of a shear force pump. Is there a way that he could illustrate this more 
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fully? One possible explanation for this phemonenon at the tip is that the 
much higher loss there would result in too low a pressure — were it not for 
the action of the blades in maintaining radial equilibrium. Accordingly, 
the outlet tangential velocity component becomes quite large at the tip in 
an effort by the machine to produce more work there to offset this loss. 
Now, regardless of the blade-to-blade and hub-to-shroud velocity profiles, 
is not this process the result of guidance by the blades? Is the tip axial 
velocity component still connected with the tangential component 
through the blade angle there? 

F. F. ANTUNES (Ingorsoll Rand Inc.) : Since the testing was done as a 
free impeller (without a stator or guide vanes) , there is some question as 
to whether the three-dimensional effects are greatly exaggerated from 
that of an inducer with stator downstream. It is suggested that a reason- 
able stator be added and velocities remeasured. 


LAKSHMINARAYANA (author) : I would like to thank the dis- 
cussors for their helpful comments and criticisms. 

Professors Wislicenus and Horlock have commented on the mysterious 
effect ; i.e. , large change in angular momentum as the flow proceeds down- 
stream. The author explained that this is due to large redistribution of 
mass flow that occurs downstream of the inducers. This can be explained 
from the following calculations. To chock the accuracy of the measure- 
ments, area averaged axial velocity V z and mass averaged Ve are com- 
puted at the measuring stations using the following equations. 


Vm 

Ut 


-J-f 1 

1-RS 4 


Vj 

u, 


R dR 


V, 

U t 


f 1 Ve V* 
L Ut Ut 


R dR 


f 1 — 

Rh U, 


R dR 


The values of V z /U t , Ve/U t at various measuring stations (see fig. 10) 
are shown in table D-I. 

As can be seen from table D-I, the accuracy of V z measurements, 
which is within 4 percent, is better than expected in such a flow. The 
change in mass averaged Ve is small compared to local changes observed 
at various stations (compare fig. 11 and table D-I). This confirms the 
author’s conclusions that local Ve changes are brought about mainly by 
flow redistribution. The mass averaged Ve shows a consistent trend in 
that it decreases continuously downstream, the decrease being small 
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Table D-I . — Values of V z /Ui, Ve/Ut at Various Stations 


Station number 

4 

5 

6 





; 

0.076 

0.0775 

0.0745 

Vo 

0.346 

0.326 

0.306 

u, 





everywhere except between stations 6 and 7. The reasons for this are 
mentioned by Wislicenus. Some of the discrepancy in T'« may also be due 
to experimental inaccuracy in measuring the flow in the backflow region, 
whose extent extends up to 20 percent of annulus at station 7. 

To check whether a simple throttle valve used in these experiments has 
any effect on inducer flow, especially in inducing backflow at the hub 
region, measurements were carried out with a properly designed throttle 
as shown in the insert in figure D-l. The tests were carried out with a 
three-bladed inducer, since this geometry was readily available. Please note 
that all the measurements reported in this paper are earth'd out in a 
four-bladed inducer and hence velocity profiles plotted in figures D-l and 16 
are for different solidity inducers. The results shown in fig. D-l indicate 
that the throttle has no ('fleet on flow separation at the blade; trailing 
edge and its influence is felt only beyond station 6, which is far down- 
stream of the inducer. Thus, it is concluded that the throttle has no effect 
on flow separation observed at stations 3, 4, and 5, even though the extent 
of backflow is somewhat reduced at station 6. 

The author agrees with Mr. Crouse with regard to the source of the 
losses. Knowledge of the source of losses is extremely important in the 
design of inducers where incorporation of viscous effects is a necessity. 
The group at Penn State is now engaged in understanding the various 
sources of losses. This would lead to an analytical model which recognizes 
various mechanisms responsible for flow departure and losses. 

In regard to Professor Horloek’s comment oirthe magnitude of the body 
force term due to viscous losses, I have proved below that its magnitude 
is extremely small and its neglect in the equations of motion is justified. 

Since the radial equilibrium equation used in this paper is passage 
averaged, it is necessary to find the average (blade to blade) of the body 
force term due to viscous losses. This term is thus equal to 
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1 d /AP„\ 

^3 tan ^(T) (3?) 

(from equation (13)). 

In deriving equation (37) , it is assumed that the streamwise gradient of 
local loss in P„ is nearly the same as the corresponding gradient of average 
loss in P 0 . Using equation (21), it can be proved that 

— 1 —dW 

Fd r ~~ X* It' — — tan a (38) 

o OX 



Figure D-l. — Radial variation of V,/Ut and \p t at various axial stations with simple 
(old) and aerodynamically designed (new) throttle (three-bladed inducer). 
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since 

L 1 

d h cos 0 

and 

d /aF 0 \ \r — dW , F 2 d\ R 
dr\pj r, dR 2 r, dR 

Hence the ratio of the two terms is 


F ir 


dW 


r,X« tan a 

dx 


d 

dr 



d W W 3X/A 
R dR + ~2dR/ 


(39) 


(40) 


This ratio is evidently very small as illustrated in this example (using 
measured values) . 


At/f = 0.8, tan a = 0.36 


X.-0.19, ^-2.1 


(from fig. 13) 
(from fig. 9) 


— dW 

W = 0.GU„ — = 0.6(7, 
dR 


(from fig. 14) 


dW -0.05(7, 
dx r t 


(fig. 9 and eq. (27) ) 


Substituting these values in equation (40) , we get 

~F~ 

— ^1.5X10~ 3 

d_ /AAA 
dr\p) 


Hence the body force term neglected in my analysis is extremely small 
compared to d/dr(AP 0 /p) and its neglect is certainly justified. Even at 
the tip, where 31 V/dR is negative, the ratio is extremely small since 
dXrt/dR is very large and a is very small. 

The other comment by Professor Horlock concerns the choice of the 
parameter </> in the loss correlation. The physical reasoning behind this 
choice is that at lower 4> (which means larger stagger angles) narrower 
blade passages result and thus the flow in these passages becomes fully 
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developed sooner and hence viscous effects become larger. I don’t agree 
with Professor Horlock that </> - ^design should be the parameter. Unlike a 
compressor aerodynamieist, an inducer hydrodynamicist is faced with the 
dilemma of evaluating the losses even at the design conditions where the 
frictional losses are several orders of magnitude higher than those in a 
compressor. It is the intention of the author to understand the design 
flow first and then attempt to correlate the flow losses at off-design 
conditions. 

The physical flow process occurring within the inducer is more com- 
plex than explained by Dr. Cooper. Larger increase in V e and absolute 
head rise is not due to guidance of the blade. The flow deviation angles 
are found to be large at every radial location. As explained in this paper, 
the radial velocities are large and are of the same order of magnitude as 
the axial components. Under these circumstances, the inducer should be 
treated as a mixed flow pump. The meridional component (in the viscid 
flow) through the inducer can be divided into axial and radial com- 
ponents. The axial component can be treated in the manner described by 
Dr. Cooper, but the radial component of the flow makes the inducer 
behave as a shear force pump. The author is presently investigating the 
possibility of predicting the flow by this approach. 

With regard to the comment by Mr. Antunes, it is very doubtful that a 
stator or guide vane would improve the flow through the inducer. On the 
other hand, large flow redistribution that occurs downstream of the 
inducer may be controlled, to some extent, by means of a stator. From 
physical and mathematical considerations, it is very hard to see how a 
stator would change the three-dimensional viscous character of the flow 
inside the inducer. 



